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ABSTRACT 

PSR J1744— 3922 is a binary pulsar exhibiting highly variable pulsed radio emission. We report 
on a statistical multi-frequency study of the pulsed radio flux variability which suggests that this 
phenomenon is extrinsic to the pulsar and possibly tied to the companion, although not strongly 
correlated with orbital phase. The pulsar has an unusual combination of characteristics compared 
to typical recycled pulsars: a long spin period (172 ms); a relatively high magnetic field strength 
(1.7 x 10 10 G); a very circular, compact orbit of 4.6 hours; and a low-mass companion (0.08 Mq). 
These spin and orbital properties are likely inconsistent with standard evolutionary models. We find 
similarities between the properties of the PSR J1744— 3922 system and those of several other known 
binary pulsar systems, motivating the identification of a new class of binary pulsars. We suggest that 
this new class could result from either: a standard accretion scenario of a magnetar or a high-magnetic 
field pulsar; common envelope evolution with a low-mass star and a neutron star, similar to what is 
expected for ultra-compact X-ray binaries; or, accretion induced collapse of a white dwarf. We also 
report the detection of a possible K' = 19.30(15) infrared counterpart at the position of the pulsar, 
which is relatively bright if the companion is a helium white dwarf at the nominal distance, and discuss 
its implications for the pulsar's companion and evolutionary history. 

Subject headings: binaries: pulsars: individual (PSR J1744— 3922) — binaries: pulsars: evolution 
(PSR J1744-3922) 



1. INTRODUCTION 

A mid-Galactic latitude pulsar survey with the Parkes 
Radio Telescope ()Crawford et alj l2006f) detected three 
new pulsars in binary systems, none of which easily 
fits within the standard evolutionary scenarios proposed 
for the majority of recycled pulsars. One of them, 
PSR J1744— 3922, was independently discovered during 
the repro cessing of the Parke s Multibeam Pulsar Sur- 
vey data ([Faulkner et al.ll200l ). This 172-ms binary pul- 
sar has a relatively high surface dipole magnetic field 
strength (B = 3.2 x 10 19 (PP) 1 / 2 G = 1.7 x 10 10 G) sug- 
gesting it is mildly recycled. However, it appears to have 
a very light companion (minimum mass 0.085 Mq) in 
a tight and nearly circular 4.6-hr orbit (see Table [T|). 
This type of orbit and companion are typical of those of 
fully recycled pulsars (which we define as pulsars with 
P < 8ms and B < 10 9 G) with He white dwarf (WD) 
companions. Thus, why PSR J1744— 3922 escaped being 
fully recycled is a puzzle. 

In addition to this atypical combination of spin and or- 
bital properties, PSR J1744— 3922 exhibits strong pulsed 
radio flux modulations, making the pulsar undetectable 
at 1400 MHz for lengths of time ranging from a few tens 
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of seconds to t e ns of minutes. It has been suggested by 
iFaulkner et alj (|2004f ) that this behavior might be the 
nulling phenomenon seen in a handful of slow, isolated 
pulsars. Nulling is a broad-band, if not total, interrup- 
tion of the radio emission for a temporary period of time. 
On the other hand, although nulling could affect pulsars 
in binary systems as well, many binary pulsars vary due 
to external effects such as eclipses. Such external effects 
might explain PSR J1744-3922's variability as well. 

In this paper, we first report on multi-frequency obser- 
vations of this pulsar that suggest the radio variability 
is not intrinsic to the pulsar. However, our analysis does 
not show strong evidence of a correlation between radio 
flux and orbital phase, as one might expect from tra- 
ditional eclipse-like variability. We then report on our 
infrared search for a counterpart of the companion using 
the Canada- France-Hawaii Telescope (CFHT). This ob- 
servation has identified a K' = 19.30(15) star at the po- 
sition determined by the radio timing observations. Fi- 
nally, we examine why other properties of this pulsar 
make it incompatible with standard evolutionary scenar- 
ios, and identify a few other systems which have similar 
characteristics. The addition of PSR J1744— 3922 as an 
extreme case among this group motivates us to identify 
a possible new class of binary pulsars. We propose sev- 
eral possible evolutionary channels which might produce 
members of this class and explain how the nature of the 
companion to PSR J1744— 3922 could be used to con- 
strain the origin of these systems. 

2. PULSED RADIO FLUX VARIABILITY 

The observed average pulsed radio emission from a pul- 
sar can fluctuate for several different reasons. These 
inclu de effects fr om the pulsar itself, as in nulling 
(e.g. lBackenfl970l ). its environment, as in eclipsing bi- 
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nary pulsars (e.g. iFruchter, Stinebring fe Tavlorjj l988'). 
or th e interstellar medium, as in scintillation (e.g. lRickettl 
I1970H . In the case of PSR J1744-3922, during a typi- 
cal 1400 MHz observation, the radio emission seems to 
turn on and off randomly on timescales varying from 
tens of seconds to tens of minutes (see Figure [T] for 
sample folded profiles). In many instances, the pulsar 
is undetectable through entire observations (~ 15min). 
In a previous analysi s of these radio flux modulations, 
iFaulkner et al.l (|2004f ) concluded on the basis of obser- 
vations at 1400 MHz that PSR J1744-3922 is probably 
a pulsar experiencing pulse nulls. As nul ling is a broad- 
band phenomenon (e.g. lBartel et alll98lf l. we decided to 
investigate the frequency dependence of the fluctuations 
in order to test the nulling hypothesis. 

2.1. Data and Procedure 

The work we report is based on an extended 
dataset c ombining both the d ata repo rted indepen- 
dently by IFaulkner et al.l ((2004) and by iRansom et alJ 
(|in prep. I ). A total of 112 radio timing observations of 
PSR J1744-3922 were made at the 64-meter Parkes Tele- 
scope and the 100-meter Green Bank Telescope between 
2003 June and 2006 January. Relevant details for the 
current study are summarized in Tableland we refer to 
the above two papers for more details about the obser- 
vational setups and timing results. 

For the purpose of studying the radio emission variabil- 
ity, we made time series of the pulsed flux intensity. We 
dedispersed the data at the pulsar's dispersion measure 
(DM) of 148.5 pc cm~ 3 and then folded the resulting time 
series in 10-s intervals u sing the timing ephemerides from 
IRansom et al.l (|in prep J ). For each observation, the pulse 
phase was determined from the profile averaged over the 
entire observation. We fit each 10-s interval of the folded 
pulse profile with a constant baseline plus a Gaussian of 
variable amplitude having a fixed width at the predeter- 
mined phase. A Gaussian FWHM=0.01964P, where P 
is the pulse period, nicely fits the profile averaged over 
many observations in the frequency range 680-4600 MHz. 
Errors on the best-fit amplitudes returned by our least- 
square minimization procedure were scaled under the as- 
sumption that the off-pulse region RMS represents the 
total system noise. Although no flux standard has been 
observed, we obtained pulsed flux density estimates by 
using the radiometer equation and equating predicted 
noise levels to the off-pulse RMS level. We also accounted 
for the offset between the telescope pointing and the real 
position of the source (since in the early observations the 
best-fit timing position had not yet been determined) 
by approximating the telescope sensitivity to be an az- 
imuthally symmetric Gaussian having a FWHM corre- 
sponding to the radio telescope beam size, which is 8.8' 
and 13.8' at 1400 MHz for GBT and Parkes, respectively. 

In this way, we generated flux time series for all 112 
observations of PSR J1744-3922 (see Figure [2] for ex- 
amples). As we describe next, these results show that 
scintillation and nulling are highly unlikely to be the ori- 
gin of the observed variability. 

2.2. Radio-frequency-dependent Variability 

As IFaulkner et alJ (|2004fl discussed previously, inter- 
stellar scintillation (ISS) is unlikely to be the source of 



fluctuations in PSR J1744-3922. Scintillation is pro- 
duced by a diffractive scattering medium along our line 
of sight and the typical diffractive scintillation t imescale 
can be expressed as (see lCordes fc Rickettlll998t) : 

At d = 2.53 x 10 4 ^=^ s, (1) 

with D the distance to the source in kpc, Av^ the decor- 
relation bandwidth in MHz, v the observed frequency in 
GHz and Viss the velocity of ISS diffraction pattern in 
kms . For the 1400MHz observation shown in Figure 
[H for example, the NE20 01 model for the Galac tic distri- 
bution of free electrons (jCordes fc Laziol 120021 ) predicts 
&v d = 0.01 MHz in the line of sight of PSR J1744-3922 
at a distance of 3.0 kpc estimated from the DM. Viss is 
typically dominated by the pulsar velocities, which are 
in the range 10-100 kms -1 for most binaries. Therefore, 
we estimate the scintillation timescale to be of the order 
of a few seconds to a minute at most. This could be com- 
patible with the fast flux variations but can certainly not 
explain the extended periods where the pulsar goes unde- 
tected. Perhaps most importantly, strong scintillations 
will be averaged away since typical observing bandwidths 
are much larger than the decorrelation bandwidth, and 
therefore contain many "scintles" . Such averaging effec- 
tively rules out the ISS hypothesis. 

Another possibility is that the flux modulation is re- 
lated to intrinsic nulling of the pulsar. Based o n observa- 
tions at 1400 MHz only. IFaulkner et all 1)20041 ) identified 
it as the most likely explanation. Only observed in old, 
isolated pulsars so far, (though nothing prevents a binary 
pulsar from being a miller) nulling is a broa d-band, if not 
total , interruption of the radio emission (e.g. lBartel et al.l 

[Mil). 

Considering the fraction of observations with no detec- 
tion of radio emission at various frequencies (see Table[3]), 
we note qualitatively that PSR J1744— 3922 is regularly 
undetectable at low frequencies but easily detectable at 
high. For instance, we obtained four Parkes observa- 
tions at 680 and 2900 MHz simultaneously in which the 
pulsar is detected twice at 2900 MHz while remaining un- 
detected at 680 MHz. In addition, seven long GBT ob- 
servations centered at 1850 and 1950 MHz show highly 
variable emission but little evidence that the pulsar ever 
disappears completely (see Figure [5]). Clearly, however, 
observations could be biased by the relative instrumental 
sensitivity in each band and by the intrinsic spectrum of 
the pulsar. 

To investigate the effect of instrumental sensitivity and 
spectral energy distribution, we analysed the pulsed flux 
densities at different frequencies. Measured values were 
estimated using the radiometer equation as explained in 
Section 12.11 and are d isplayed in Table |U We note that 
IFaulkner et alJ (|2004h reported a different flux density 
than ours at 1400 MHz (0.20(3) vs. 0.11(3) mJy, respec- 
tively). The discrepancy could be because the average 
flux density changes depending on the amount of time 
PSR J1744-3922 spends in its "bright state" during an 
observation. The large standard deviation (0.16 mJy, see 
TableH]) at this frequency suggests that by restricting the 
calculation to observations for which PSR J1744— 3922 is 
nicely detected, a higher flux value can be obtained. 

For the simultaneous observations at 680 and 
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2900 MHz in which the pulsar was not detected at 
680 MHz (see Table [2]), we can put an interesting ap- 
proximate lower limit on the spectral index, a, of the 
pulsar if we assume that the minimum detectable pulsed 
flux at 680 MHz is an upper limit to the pulsed flux at 
this frequency: 

log(S , 2900/S , 680) n ro s 

&>-, — -. -, -~0.17. (2) 

log (^2900/^680 J 

Such a value is extremely flat compared to that of 
the average populatio n of pulsars, which has a spectral 
index of -1.8 ± 0.2 (|Maron et al.l 12000). Thus either 
PSR J1744— 3922 has a spectrum very different from 
those of most pulsars and/or the flux variability is in- 
trinsically frequency dependent. 

The distribution of pulsed flux density values for all 
1400 and 1950 MHz observations is shown in Figure[31 It 
is clear from the distribution at 1400 MHz that the nu- 
merous non-detections are responsible for the peak be- 
low the sensitivity threshold. We also observe that the 
1400 MHz flux density has a higher average value and is 
much more variable than at 1950 MHz. This frequency 
dependence suggests some sort of scattering mechanism 
with the unscattered flux level much higher than the ob- 
served average flux at 1400 MHz, and argues against it 
being classical nulling. Therefore, we conclude that this 
unknown mechanism affecting the lower frequency flux 
might be responsible for the apparent flat spectrum de- 
rived from the simultaneous 680-2900 MHz observation. 

Recent simultaneous multi-frequency observations of 
PSR B1133+16, a well-established nuller, show that sin- 
gle pulse nulls do not always hap pen simultaneously at 
all frequencies (|Bhat et al.1 l2006f) . They also observe, 
however, that the overall null fraction does not present 
any evidence of frequency dependence, which might mean 
that sometimes nulls are simply delayed at some frequen- 
cies. In the case of PSR J1744-3922, the S/N limits 
us to consider the pulsed flux over times correspond- 
ing to many pulse periods only. Therefore, the kind 
of non-simultaneo us, frequency-dependent effect seen by 
iBhat et "ail (|2006f ) is not relevant to our analysis and thus 
we expect the variability to be independent of frequency 
if really caused by nulling. 

Although our flux measurements at other frequencies 
are not simultaneous, we can assume they are good sta- 
tistical estimates of the normal flux of the pulsar and 
use them to characterize its spectrum. In an attempt to 
estimate an unbiased spectral index, we can use the ap- 
proximate maximum flux value at each frequency. The 
1400, 1850, 1950, 2900 and 4600 MHz data give a spec- 
tral index lying between —1.5 and —3.0, which is similar 
to many known pulsars. However, it seems that the flux 
at 820 MHz is much smaller than expected from a sin- 
gle power-law spectrum. Since flux variations are very 
important at low frequency and we only have a single 
detection at 820 MHz, the reported maximum value is 
probably not representative of the real flux of the pulsar 
at this frequency. 

In summary, the facts that the pulsar radio emission 
rarely drops below our detection threshold at 1950 MHz 
and that the radio variability is frequency dependent, 
demonstrate that the flux modulation is probably not 
classical nulling. A non-nulling origin for the fluctuations 



at 1400 MHz also explains why PSR J1744-3922 does 
not fit in with expectations based on the correl ations 
obser ved between null fraction and spin period (Biggs 
1992), and between nul l fraction and equivalent pulse 
width (jLi fc Wanglll995h . In comparison with millers, it 
has one of the smallest spin periods and a small pulse 
width (~ 3.4ms), but one of the largest "null" fractions 
(~ 60% at 1400MHz). Since this pulsar is in a tight bi- 
nary system, the possibility of influence by its companion 
is therefore an important alternative to consider. 

2.3. Orbital Correlation Analysis 

Even though a quick examination of the time se- 
ries confirms that the flux decreases observed for 
PSR J1744— 3922 are not due to systematic eclipses of 
the pulsar by its companion, a more subtle orbital cor- 
relation could exist. To search for such an effect, we ask 
whether or not the pulsar is more likely to be detected 
at a particular orbital phase. For this analysis, we folded 
the time series in 1-min intervals, and defined the pulse 
as detected if the best-fit Gaussian amplitude was greater 
than its la uncertainty. In order to limit spectral effects, 
we restrict the analysis to observations covering the range 
1237.5-1516.5 MHz at Parkes 7 and 1404.5-1497.5 MHz at 
GBT. 

The results of this analysis are shown in Figure [4^. 
The histogram represents the fraction of detected pulses 
with respect to the total number observed, as a function 
of orbital phase. Errors were determined using Poisson 
statistics, implicitly assuming our assigning of each inter- 
val as a detection or non-detection is accurate. There is 
a suggestion that PSR J1744— 3922 is more often unde- 
tectable between phases ~ 0.05 — 0.45. The best-fit con- 
stant model gives a x 2 /9 = 5.85 (the histogram has 10 
orbital phase bins), which, if correct, would be highly sig- 
nificant. To test the accuracy of our errors, we performed 
a Monte Carlo simulation consisting of 10000 trials where 
we assigned each measurement a random orbital phase. 
The mean x 2 /9 = 0-37 with a standard deviation of 0.18, 
suggesting our error estimates are significantly overesti- 
mated. 

To investigate further, we performed an analysis sim- 
ilar to the previous one but for the pulsed flux density 
measured at each orbital phase averaged over all obser- 
vations. We selected two subsets of data: the Parkes and 
GBT observations made at 1400 MHz and the GBT 1950 
MHz observations. The first one includes 101 observa- 
tions which are on average ^15 minutes long, but there 
are a few observations which are significantly longer. The 
latter subset includes four observations, two of which 
have full orbital coverage, one covering ~ 75% of the 
orbit and the last one ~ 40%. Observations with no 
detection of PSR J1744— 3922 are assigned upper limit 
values of three times the background noise level (which 
is very small compared with the average pulse of the 
pulsar when it is on). Errors in each bin of the his- 
togram are estimated from the RMS of the individual 
values in each orbital bin. Results are plotted in Figure 
Qi). For the 1400 MHz data, a fit to a constant line has 
a x 2 /9 = 15.58, and for the GBT 1950 MHz data, we 
find a % 2 /9 = 5.59. Randomizing the individual data 

7 This is the common range of the observing modes centered at 
1375 and 1400 MHz. 
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points and folding them resulted in a x 2 /9 — 1-02 at 
both frequencies, suggesting our error estimates are rea- 
sonable. Although this analysis strongly rules out the 
constant model for our folded light curves, the shapes 
of these curves at 1400 and 1950 MHz are not consistent 
with each other. 

This could be a result of the paucity of observations 
(typically 3-5) at any given orbital phase. Therefore, 
random fluctuations in the flux on timescales of tens of 
minutes (which we see in the time series) would likely not 
be averaged out. To test whether or not the significant 
deviation from a constant model depends on the par- 
ticular phasing of our observations, we again performed 
Monte Carlo simulations in which 10000 trial histograms 
similar to those shown in Figure [4] were generated from 
the real data by adding a random orbital phase shift at 
the beginning of each time series and then determining 
the x 2 /9 value for a constant model. From the result- 
ing distribution of x 2 /9 values, we estimate the chance 
probability of obtaining the particular x 2 /9 values ob- 
tained, or higher, using the real orbital phases. For 
the on-off analysis, the probability is 0.102, or a formal 
x 2 /9 = 1.63. For the 1400 and 1950 MHz flux density 
light curves, the probabilities are 0.015 (x 2 /9 = 2.28) 
and 0.212 (x 2 /9 = 1-34). This suggests there may be 
some correlation with the orbit, but the shape of our 
folded light curves are still dominated by more stochas- 
tic flux variations given our limited data set. We would 
expect a standard eclipse mechanism to make the pulsar 
dimmer when the companion is in front of it at phase 
0.25, which does not seem to be the case. There may 
be large orbit to orbit variations in the eclipse depths, 
durations, and phases. This kind of behavior has been 
observed in other binary pulsar sys tems such as Ter5A , 
Ter5P, Ter5ad and NGC6397A (seelRansom et al.l l2005l 
iHessels et al.l l2006t iD'Amico et al.ll2001l for examples), 
but a much larger data set would be required to obtain a 
reliable average light curve to show if this is the case for 
PSR J1744-3922. It is still possible that the short time 
scale pulsed flux variations tend to group in "events" dur- 
ing which the pulse gets dimmer. These "events" may 
last for a significant fraction of the orbital phase caus- 
ing the apparent marginal orbital correlation given our 
limited statistics. 

2.4. Accretion and mass loss limits 

Many other systems are known to exhibit strong 
flux radio variations for which the orbital dependence 
is well establishe d. One of them is PSR B1718-19 
(|Lvne et al.lll993l ). Interestingly, this pulsar has a low- 
mass companion and orbital properties similar to those 
of PSR J1744-3922 and is also harder to detect at low 
frequency. At 408 and 606 MHz, PSR B1718-19 gets 
so dim that it is barely detectable in spite of good sen- 
sitivity during a large part of the orbit. On the other 
hand, at 1404 and 1660 MHz the orbital modulation of 
the average flux is much less important. This flicker- 
ing is probably made by material left over by the wind 
of the companion, a bloated m ain sequence (MS) star 
(jJanssen fc van Kerkwiikll2005f ). Although the compan- 
ion of PSR B1718— 19 is not large enough to come near 
to filling its Roche-lobe, this could happen in a tighter 
binary system like that of PSR J1744-3922. In this 
case, some kind of tidal stripping could be occuring, 



leaving material around the system. This could explain 
why the pulsar does not disappear at conjunction like 
PSR B1718— 19 does, but in a more stochastic way. 

If the companion is losing mass, one might expect to 
observe small ch anges in the orbital par ameters. From 
radio timing (see lRansom et aLllin prepj ). we can set an 
upper limit of \P or b\ < 2 x 10 _10 ss~ 1 and |i| < 7 x 
10 -12 lt-ss -1 on the rate of change of the orbital period 
and the rate of change of the projected semi-major axis, 
respectively. We can use the latter quantity to evaluate 
the implied mass loss limit. For a circular orbit, which 
is a good approximation here, we can express th e rate of 
change of the semi-major axis as (IVerbuntlll993l) : 



a J M c V M„ 



(1-/?)M C M 
a(l - (3)- 



2(M P + M c ) 



M p - 

where x = a sin i, x — a sin i , M c and M p are the mass of 
the companion and the pulsar, respectively, J is the total 
angular momentum of the system, (3 is the fraction of 
mass accreted by the pulsar and a is the specific angular 
momentum of the mass lost in units of the companion 
star's specific angular momentum. 

For the case in which the total orbital angular mo- 
mentum of the system is preserved (J = 0) we can see 
that mass loss from the companion 8 (M c < 0) would 
necessarily lead to a widening (a /a > 0) of the orbit 
if: 1) the mass transfer is conservative ((3 — 1), or 
2) the mass transfer is non-conserva tive (f3 < 1) and 
a < 1 + M C /(2M P ) fsee lVerbuntl[l993L for more details). 

By considering the conservative case, in which |M C | = 
\M p \, we obtain an upper limit on a possible mass ac- 
cretion rate by the pulsar \M p \ < 3 x 10~ 12 M Q yr _1 . 
This would be even lower for the non-conservative case. 
Accretion onto the pulsar is possible if the corotation 
radius: 

r co ~ BOOP^M^ 3 km, (4) 
is larger than the magnetospheric radius of the pulsar: 



^mag 



800 



R4 pl2 

Mi. 4 M£ 12 



1/7 



km, 



(5) 



where P psr and B psr are, respectively, the spin pe- 
riod and the surface dipole magnetic field in units of 
PSR J1744-3922's P and B (172 ms and 1.68 x 10 10 G); 
Mi. 4 is the mass of the pulsar in units of 1.4 M Q ; Rio is 
the radius of the pulsar in units of 10km; and Af_i2 is 
the accretion rate in units of 10 -12 M Q yr _1 . Given the 
upper limit on the mass accretion rate, and the likely 
conservative Rio — 1 value for the neutron star radius, 
we find that r rnag (> 600 km) is likely larger than r co 
(~ 500 km) . This argues against any significant accre- 
tion occuring in the system. 

That no significant accretion is occuring is also sup- 
ported by an XMM-N ewton observation that allows 
iRansom et al.1 ( in prep. I ) to put a conservative upper 
limit ~ 2 x 10 ergs" 1 on the unabsorbed X-ray flux 
from accretion in the 0.1-10keV range. Assuming 

8 Here we implicitly assume that M c < M p . 
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with a conversion efficiency of accretion energy into ob- 
served X-rays r\ = 0.1, we get M p < 2 x 10~ 14 M yr" 1 
for accretion at the surface of the pulsar and M p < 
2.4 x 10~ 12 Mq yr -1 if accretion is limited to the bound- 
ary of the magnetospheric radius. Therefore, it appears 
unlikely that PSR J1744-3922 is accreting and, if the 
companion is losing mass, it is probably expelled away 
from the system. 

3. INFRARED OBSERVATIONS 

The radio flux variability, which might be due to mate- 
rial leaving the surface of the companion, and the atypi- 
cal evolutionary characteristics of the pulsar (see Section 
13]) can be investigated further by observing its companion 
at optical or near- infrared wavelengths. We imaged the 
field of PSR J1744-3922 at K'-band on the night of 2005 
April 19 using the Canada-France-Hawaii 3.6-m Tele- 
scope (CFHT) at Mauna Kea. The telescop e is equipped 
with the Adaptive Optic Bonette (AOB, iRigaut et al.l 
1998), which provides good correction for atmospheric 
seeing, and KIR, the 1024x1024 pixel HAWAII infrared 
detector with 0V0348 pixel scale. The total integration 
time was 30sx59 integrations = 1770s. 

We substracted a dark frame from each image, and 
then constructed a flat-field image from the median of 
the science frames. The images were then flat-fielded, 
registered and stacked to make the final image. The final 
stellar profile has a FWHM of 0'.'17, degraded somewhat 
from the optimal correction provided by the AOB system 
due to a high airmass (~ 2.0) and poor natural seeing. 
Figure [5] shows the final image we obtained. 

We analyzed th e final CFHT image using the standard 
routine daophot (|Stetsonlll987h for PSF fitting photom- 
etry, and calibrated the image using the standard star 
FS34 (|Casali fc Hawardedfl992l) . To find the photomet- 
ric zero point, we performed photometry on the standard 
star with a large aperture containing most of the flux, 
and applied an aperture correction for the PSF stars in 
the science image. Careful calibration of measurement 
errors has been done by adding artificial stars of known 
magnitude to blank parts of the image and then mea- 
suring their magnitude through the PSF fitting process 
along with the real stars. Thus, errors on the magnitude 
returned by daophot can be rescaled by calculating the 
standard deviation for the added stars. 

We found the astrometric solution for the image by 
cross-identifying five stars with the 2MASS catalogue 
(|Skrutskie et aLlfeOOGj ). fitting for scale, rotation and dis- 
placement. The final astrometric uncertainty is 0'.'34 
at the 3<7 confidence level. This value depends on the 
matching to the reference stars because the error on 
the radio timing position of PSR J1744— 3922 is negli- 
gible, (~ 0'.'03). The final image is displayed in Figure 
with the positional error circle centered at the radio 
position of PSR J1744-3922: a = 17 h 44 m 02!667(l) and 
S = -39°22'21'/52(5). Only one star falls inside this cir- 
cle, and for this we measure K' = 19.30(15). We observe 
that, above the 3cr detection limit, the average stellar 
density is 0.079 arcsec -2 and hence the probability of a 
star falling in the error circle is only 2.9%. Due to its 



positional coincidence and the low probability of chance 
superposition, we henceforth refer to this object as the 
possible counterpart to PSR J1744-3922. 

Unfortunately our near-infrared observation does not 
tightly constrain the nature of the companion to 
PSR J1744— 3922, mainly because of the uncertainties in 
the distance to the system and in the companion mass, as 
well as the fact t hat its temperature is u nknown. Assum- 
ing the NE2001 (jCordes fc Lazioll2002f ) electron density 
model is correct, we infer a DM distance of 3.0 ± 0.6 kpc 
and hence a distance modulus ranging from 11.9 to 12.8. 
U sing the three-dimens ional Galactic extinction model 
of iDrimmel et al.l {2003) we get a value of Ay — 1.9 for 
a distance of 3.0 kpc. Converting 9 t he inferred extinc- 
tion to K' band gives A K , ~ 0.21 fsee lRieke fc Lebofskvl 
[19851 for conversion factors). Therefore the estimated 
absolute magnitude of the counterpart lies in the range 
M K > a [6.1, 7.4]. 

We can evaluate how probable it is that the compan- 
ion is a He WD, a typical low-mass companion in bi- 
nary pulsar systems, since WD cooling models can put 
restrictions on the stellar mass and cooling age, given 
an observed flux. Figure [6] shows the absolute K' mag- 
nitude as a function of cooling age for He WDs of dif- 
ferent masses. These cooling tracks were made by us- 
in g WD atmosphere mo dels based on the calculations 
of iBereeron et all (119951) , and thereafter improved by 
iBergeron et al.l (|2001h and iBergeronl (|2001l ) in combina- 
tion w ith evolution sequences calculated bv lDriebe et alJ 
(1999). Although the mass range of models available to 
us does not go below 0.179 M Q , we can deduce that to 
be so luminous, a He WD would need to have a very low 
mass and a cooling age much lower than the character- 
istic age of the pulsar (1.7 Gyr). Even if the mass were 
equal to the lower limit of 0.08 M Q derived from radio 
timing, it seems unlikely that such a companion could 
be as old as the pulsar characteristic age. Therefore, if 
the companion is indeed a He WD, the pulsar's spin pe- 
riod must be close to the equilibrium spin frequency it 
reached at the end of the recycling process. In this case, 
its characteristic age is an overestimate of its true age. 

Another possibility is that the companion is not a He 
WD. If the companion is a low-mass main sequence (MS) 
star, then the minimum mass required to match the lower 
limit on the absolute K' flux is ~ 0.25 M Q , regardless 
of the pulsar age. Such a companion mass requires a 
favorable face-on orbit but this cannot be ruled out from 
the near-infrared and radio data. On the other hand, a 
lower- mass bloated MS star could be equally as bright, 
so is also a possibility. 

We cannot constrain the nature of the counterpart very 
well from a measurement in a single near-infrared filter. 
Ideally, obtaining a spectrum could yield: 1) a precise 
determination of the nature of the counterpart, 2) orbital 
Doppler shift measurements of the spectral lines which 
can prove the association as well as determine the mass 
ratio, 3) in the case of a white dwarf, an estimate of its 
cooling age from spectral line fitting. 

4. BINARY PULSAR EVOLUTION 



9 For simplicity, and because the error on the magnitude is dom- 
inated by the distance estimate, we hereafter assume that K and 
K' are similar. 
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The sporadic radio emission from PSR J1744— 3922 is 
not the only indication that there is something unusual 
about the pulsar's interaction with its companion. The 
pulsar is in a tight and low eccentricity (e < 0.001) or- 
bit with an apparently very light companion having a 
minimum mass of 0.08 Mq. On the other hand it has a 
relatively large surface magnetic field (1.7 X 10 10 G) and 
an extremely long spin period (172 ms) compared with 
othe r binary pu l sars having low-mass com panions (see, 
e.g.. IStairsll2004t Ivan Kerkwiik et al.ll2004 for recent re- 
views). These properties, along with the relatively bright 
near-infrared counterpart of the companion, make it un- 
usual and suggest it evolved differently than most binary 
pulsars. 

In binary pulsar evolution, the nature of the companion 
plays a key role in determining the final spin and orbital 
characteristics of the system. Most of the binary pulsar 
population consists of pulsars with low-mass companions 
in low-ecce ntricity orbit s. Following the nomenclature 
adopted bv lStairsl (|2004h . the case A evolutionary chan- 
nel results in pulsars having He WD companions. They 
are neutron stars (NS) which were spun up to very short 
periods (P < 8 ms) after accreting matter from a low- 
mass star dur ing a long and steady R oche-lobe overflow 
(RLO) phase (|Tauris fc Savoniielll999h . For reasons that 
are not yet understood, this process is responsible for the 
partial suppression of the surface magnetic field to val- 
ues of the order of 10 8-9 G. The most robust predictions 
of this model are the correlation s linking the orbital pe- 
riod to the mass of the He WD (IR appaport et al.lll995| ) 

see Figure [3 and the eccentricity to the orbital period 

Phmnevl[l992T) . 

On the other hand, the case B channel is made of 
pulsars having more massive CO WD or ONeMg WD 
companions (M c > 0.45M©). Their intermediate mass 
progenitors did not sustain a stable RLO phase, instead 
evolving in a short-duration, non-conservative, common 
envelope (CE) phase during which the pulsar spiraled 
into its companion's envelope. This process only partly 
recycled the pulsar, leading to intermediate spin peri- 
ods (P > 8 ms) and leaving a higher magnetic field 
(~10 9 - 10 G). 

In TableJH the expected properties of systems resulting 
from these two evolutionary channels arc compared with 
those of PSR J1744— 3922. Both scenarios fail to account 
for all the observed characteristics; this suggests a spe- 
cial evolution for PSR J1744— 3922. This can also be seen 
from a P—B diagram (Figure^) made for binary pulsars 
in the Galactic field having circular orbits. As opposed to 
isolated and other kinds of non-recycled binary pulsars, 
there exists a very strong relationship linking P and B 
which is presumably related to the recycling process. To 
our knowledge, such a correlation has not been reported 
in the recent literature although w as indirectly found by 
Ivan den Heuvel fe~ Bitzaraki (1995|) who reported a pos- 
sible correlation between P — P or h and B — P rb for bi- 
nary pulsars in circular orbits. The many binary pul- 
sars discovered in recent years may be making it easier 
to appreciate. In Figure [5£i we see that pulsars having 
light companions (case A) generally gather in the region 
of low magnetic field and short spin period whereas the 
case B pulsars lie in higher- valued regions. Surprisingly, 
of the six highest magnetic field pulsars, five of them, 
including PSR J1744-3922, have light companions. The 



remaining one, PSR B0655+64, is an extreme system 
associated with the c ase B channel since it has a mas- 
sive WD companion (jvan Kerkwiik et alj 120041 ) . How- 
ever, the case B channel cannot accommodate the other 
five pulsars because, assuming random orbital inclina- 
tions, a simple statistical estimate gives less than a 0.1 % 
probability for all of them to be more massive than the 
required 0.45 Mq. For PSR J1744-3922 in particular, 
the orbital inclination would need to be less than 12.5°, 
which represents a 2.5% probability. 

We also report in Table [5] the principal characteristics 
of binary pulsars that appear to be partly recycled (e.g. 
P > 8 ms and e < 0.01) but have companions likely not 
massive enough to be explained by the standard case B 
scenario. These pulsars have related properties and could 
have experienced similar evolutionary histories. Apart 
from their strange position in the P ~ B diagram, they 
also stand out when we look at the P D rb — P relation- 
ship (Figure [8Jd). In this plot, we see that pulsars having 
low- mass companions (case A) occupy the bottom re- 
gion, below P < 8 ms, and their spin periods are more or 
less independent of the orbital period. This arises from 
the fact that recycling probably sat urates for a given ac- 
cretio n rate and/or accretion mass (|Konar fe Choudhuril 
2004). On the other hand, fewer constraints exist in 
this parameter space for pulsars with massive WD com- 
panions (case B). Their short CE evolution would limit 
the recycling efficiency and thus the final parameters are 
more sensitive to the initial conditions. Finally, we ob- 
serve a third category made of relatively slow pulsars 
with low-mass companions in compact orbits. Neither 
the case A nor the case B scenarios is able to explain 
such properties, especially for the most extreme systems 
like PSR J1744-3922 and PSR B1831-00. Therefore, 
we suggest a new "class" of binary pulsars having the 
following properties: 1) long spin periods (in comparison 
to MSPs), 2) large surface magnetic fields (~ 10 10 ~ n G), 
3) low-mass companions, likely 0.08 — 0.3 M©, having 
nature yet to be determined, 4) low eccentricities, and 
possibly 5) short orbital periods (< 5d). On this last 
point, very wide orbit systems like PSR B1800— 27 and 
PSR J0407+1607 might be explained by the standard 
case A scenario in which it is difficult to achieve an ex- 
tended period of mass transfer from the companion to 
the pulsar. 

5. DISCUSSION 

The existence of another "class" of binary pulsars 
is supported by the fact that several pulsars now oc- 
cupy a region of the parameter space delimited by the 
spin period, orbital period, magnetic field and compan- 
ion mass that seems inaccessible to the standard evo- 
lutionary channels. PSRs J1744-3922, J1232-6501, 
B1718— 19 and B1831— 00 are certainly the most no- 
ticeable candidates. Although other papers also iden- 
tified that some of these pulsars have strange char- 
acteristics (see, e.g . . Ivan den Heuvel fc Bitzarakil 119951 : 
ISutantvo fc Lj||2000t lEdwards fc Bailesll2001h . there is no 
consensus on their evolutionary histories. For instance 
the case of PSR B1718— 19 was considered somewhat 
unique because it is in a globular cluster and hence has 
possibly been partly pe rturbed, or even gr eatly changed, 
by stellar interactions (lErgma et al.l 1 1996h ■ The discov- 
ery of PSR J1744-3922, in the Galactic field, is impor- 
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tant because it strengthens the possible connection be- 
tween PSR B1718— 19 and other similar binary pulsars in 
the field. Unless they do not have WD companions, these 
pulsars were at least partially recycled because tidal cir- 
cularization is needed to explain eccentricities of 0.01 
and smaller. However some of them have eccentricities 
10~ 3 — 10~ 2 , relatively large given their very short orbital 
periods, which is in contr ast to tight cas e A systems hav- 
ing smaller eccentricity (|Phmnevlll992t ). In this section 
we speculate and put constraints on some scenarios that 
may explain this possible new class of pulsars. 

5.1. Recycled High Magnetic Field Pulsar Channel 

One possibility is that pulsars like PSR J1744-3922 
may initially have had a magnetar-strength magnetic 
field (B ~ 10 14 - 10 15 G). Such a pulsar could expe- 
rience standard case A evolution involving conservative 
RLO but since the initial magnetic field is higher by 1-2 
orders of magnitude, at the end of the recycling process, 
the field might be 10 10-11 G instead, consistent with ob- 
servations. Since the recycling mechanism appears to 
strongly correlate the final magnetic field with the final 
spin period, the pulsar would also have an unusually long 
spin period as well. 

Whether this extrapolation to magnetars and high 
magnetic field pulsars is valid depends on accretion- 
induced decay models of the magnetic field strength. 
Among the proposed mechanisms is magnetic field burial, 
in which material accretes through the polar cap and, 
while piling up at the poles, exerts a latitudinal pressure 
gradient by trying to spread toward the equator. This 
effect tends to drag the field lines away from the poles, 
increasing the polar cap radius and decreasing th e mag - 
netic moment of the pulsar. iPavne &: Melatosl (|2004 ) 
show that the magnetic field would naturally "freeze" to 
a minimum stable strength once the amount of accreted 
mass exceeds some critical value, if the magnetospheric 
radius is c omparable to the size of the neutron star. Ac- 
cording to iPavnel (|2005f ). this critical mass could reach 
up to 1 M Q for a 10 15 G NS and hence magnetic field 
suppression would be difficult to achieve due to the large 
amount of accretion mass required. However, account- 
ing for other effects like the natural decay of magnetic 
field due to X-ray emission and high-energy bursts (see 
iWoods fc~T hompson 20061 for a review), might make a 
partially suppressed final magnetic field of 10 10 G plau- 
sible. 

A recycled high magnetic field pulsar is expected to 
leave a He WD companion and follow the companion 
mass-orbital period relationship as for the normal case 
A systems. PSR B1718— 19 is excluded from this kind of 
evolution because its companion is a bloated MS star. All 
the other pulsars listed in Tableware potential members, 
albeit PSRs B1800-27 and J0407+1607 would need rela- 
tively face-on orbits (i < 30 and 24°, implying 13 and 9% 
probability for randomly oriented orbits, respectively) to 
match the P or b — M c relationship (Fig. [7]). Another in- 
teresting prediction of this scenario is that, because the 
recycling process can leave pulsars with relatively long 
spin periods, they might not be very different from the 
spin periods when accretion ceased. Thus, the real age 
could be much lower than the timing-based character- 
istic age, and the WD companions would have younger 
cooling ages as well. 



The fraction of Galactic field binary pulsars in this 
class (~7/60) might naively be thought to be similar to 
the fraction of observed magnetars with respect to the or- 
dinary pulsars (~10/1500). However, the observed pop- 
ulation of magnetars suffers from severe selection effects 
because many of them appear to lie dormant, becoming 
observable for only b rief intervals, like XTE J1810— 197 
(jlbrahim e t al. 2004). Therefore, we can relate these pa- 
rameters as follows: 

fficlass _ -^obs. mag. 1 
-^binary -^radio /quies 

where N c \ ass is the number of pulsars in the new 
class, iVbj nar y is the number of binary pulsar systems, 
iV t) S mag is the number of observed magnetars, N la ^ Q 
is the number of radio pulsars and /q U i es the fraction of 
magnetars in quiescence. 

Hence we estimate that, due to quiescence, the frac- 
tion of observed magnetars with respect to the total 
population is about j^q/^q ~ 0.06. This would make 
a total population of ~ 175 magnetars which is consis- 
te nt with the possible ~ 100 in our Galaxy estimated 
bv IWoods fc Thompson! (|2006h . who used / quies =0.1. 
Our crude calculation has many caveats: Is the binary 
magnetar population similar to the binary radio pulsar 
population? Can the short lifetime of magnetars and 
high magnetic field neutron stars limit the number of 
such recycled systems? Clearly, this latter point depends 
strongly on the time evolution of the self-induced decay 
of the magnetic field, which seems to op erate on a ti me 
scale of a few tens of thousands of years (|Kaspill2004D . 

5.2. UCXB Evolutionary Track Channel 

There exists a class of neutron stars and CO/ONeMg 
WDs that accrete from light He WD donors in ultra 
tight (few tens of minutes) orbits. These X-ray emit- 
ters, known as ultracompact X-ray binaries (UCXBs), 
may be the result of wider ~ 0.5-day systems that have 
decayed due to gravitational wave radiation. A CE sce- 
nario has been proposed as a viable channel for forming 
CO/ONeMg WD - He W D and NS - He WD systems 
(|Belczvnski fc T aam 2004) . In fact such a channel would 
be very similar to the case B scenario but for an initially 
much lighter companion. For companions not massive 
enough to experience the He flash, a CE phase is possi- 
ble if the onset of mass transfer occurs late enough in the 
evolution so that the companion has reached the asymp- 
totic giant branch. In this case, the RLO becomes un- 
stable and it bifurcates from the standard case A track 
to the CE phase. 

Pulsars experiencing this evolution would be partially 
recycled, like case B systems, but they would have He 
WD companions. After this stage, only sufficiently tight 
systems with orbital periods of about one hour or less 
can evolve to become UCXBs within a Hubble time be- 
cause gravitational decay is negligible for wider orbits. If 
PSR J1744-3922-like pulsars belong to the long orbital 
separation high-end of the UCXB formation channel, we 
might expect to see more such pulsars at similar and 
smaller orbital separations than PSR J1744— 3922. It is 
possible, however, that the observed sample is biased: 
as the orbital separation decreases, wind and mass loss 
by the companion become more important and would 
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make them more difficult to detect in classical pulsar sur- 
veys conducted at low frequency (i.e. ~400 MHz) where 
eclipses are more frequent and radio emission might sim- 
ply turn o ff. The ongoing AL FA survey at Arecibo, at 
1400 MHz (|Cordes et al.ll2006[ ). could therefore find sev- 
eral new pulsars like PSR J1744-3922. Additionally, 
larger orbital accelerations make them more difficult to 
find and, since gravitational wave radiation varies as the 
fourth power of the orbital separation, their lifetimes are 
dramatically shorter. 

In such a scenario, these pulsars might be born at 
a spin period that is comparable with those of case 
B pulsars having CO WD companions, assuming that 
the short, high-accretion rate recycling would efficiently 
screen the magnetic field during the mass transfer pro- 
cess. Afterwards, because of their relatively larger mag- 
netic field, they would spin down more rapidly than case 
B pulsars. Thus, they would necessarily have longer spin 
periods and the true age is more likely to be in agreement 
with the measured characteristic age. 

5.3. AIC Channel 

Finally, a third scenario to explain the unusual prop- 
erties of PSR J1744— 3922 is the accretion-induced col- 
lapse (AIC) of a massive ONeMg WD into a NS. Most 
likely AIC progenitors would be massive ONeMg WDs 
(> 1.15 Mq) accreting from MS donors in CV-like sys- 
tems, or from red giants or He WDs in U CXB sys- 
tems (se e iTaaml 120041; H vanova fc Taamll200l for more 
details). iNomoto fc Kondol (|199lF have shown that for 
accretion rates > 0.001 Mem and/or metal-rich accreted 
mass, the ONeMg WD would collapse to a NS rather 
than explode in a supernova. More recent calcula- 
tions inc luding Coulomb corre c tions to the equation of 
state by iBravo fc Garcia-Senzl (|1999T ) demonstrate that 
AIC is possible for critical densities of the accreting 
WD core that are 30% lower than previously found by 
INomoto fc Kondol (|1991h . thus facilitating the formation 
of neutron stars through this channel. 

The properties predicted by the AIC scenario nicely 
agree with what we observe for the class we are propos- 
ing: the mass transfer prior to the NS formation would 
explain the low mass of the companion and the collapse 
is expected to be a quiet event during which almost no 
mass is lost in the system and only ~ 0.2 M is con- 
verted in binding energy into the NS. This would keep 
the final orbital period close to what it was prior to the 
collapse (which is small in most scenarios leading to AIC) 
and allow the eccentricity to be very small or, at least, 
circularize rapidly. The survival rate of such systems 
is probably higher than for standard systems which are 
more likely to be disrupted if a large amount of mass is 
lost during the supernova process. Although the initial 
properties of pulsars formed by AIC are not known, we 
may presume they resemble those of "normal" pulsars 
with magnetic fields in the 10 11-12 G range. Another in- 
teresting point to consider is that the AIC channel does 
not require a degenerate companion. As opposed to the 
other two proposed scenarios (Sections 15.11 and l5.2|) . the 
formation of a pulsar through AIC can interrupt the mass 
transfer, thus postponing further evolution. If the "evo- 
lutionary quiescence" is long enough, we would expect 
to find some non-degenerate companions around young 
AIC pulsars having circular orbits but with spin peri- 



ods and magnetic fields that are more typical of isolated 
pulsars. Finally, as the companion continues to evolve 
to fill its now larger Roche lobe, a shortened accretion 
phase might occur, thus transfering a very small amount 
of mass to the pulsar and leaving it partly recycled de- 
spite its low-mass companion. 

Some binary pulsars among the group we high- 
lighted, like PSRs B1831-00 and B1718-19, have 

been proposed as AIC candidates i n the past 

(jvan den Heuvel fc Bitzarakilll995t lErgmalll993[ ). How- 
ever, there is no firm evidence to support this. For in- 
stance, PSR B1718— 19 is pres umably a member of the 
globul ar cluster NGC 6342. As lJanssen fc van Kerkwijkl 
(2005) show, although observations suggest AIC is pos- 
sible, an encounter and tidal capture scenario can- 
not be ruled out and is very reasonable given the 
plausible globular cluster association. On the other 
hand, PSR B1718— 19 shares similarities with other 
binary pulsars in the Galactic field, especially with 
PSR J1744-3922. PSR B1718-19's younger age, larger 
magnetic field and spin period, as well as the fact that it 
has a non-degenerate companion, are all compatible with 
it being an AIC pulsar in the intermediate "quiescent" 
phase. 

In this context, the other pulsars of our putative 
class would have reached the final evolutionary stage 
and, hence, display mildly recycled properties. If so, 
PSR J1744-3922 likely has a very light He WD compan- 
ion. This might explai n the lack of tra ditional eclipses 
as in PSR B1718-19 dLvne et al.lll993h. large DM vari- 
ations as in NGC 6397A (jD'Amico et al.l [200lh . and or- 
bital pmod_d_CTiwtiyes_as_jn_^^ binary pul- 
sars (N ice. Arzoumanian fc Thorsettll2000[ ) since tidal ef- 
fects are important for non-degenerate companions. If 
the residual recycling phase left PSR J1744— 3922 with 
a long spin period, the cooling age of its hypothetical 
WD companion could be smaller than the characteris- 
tic age of the pulsar for the same reason described in 
the recycled high magnetic field pulsar scenario (see Sec- 
tion [OJ. Also, pulsars forming through AIC can have 
lower masses than those made in standard supernovae 
and since they only accrete a small amount of mass from 
their companion afterwards, they might be less massive 
than conventionally formed MSPs. 

6. CONCLUSIONS 

This paper highlights the unusual nature of the bi- 
nary pulsar system PSR J1744— 3922. The puzzling ra- 
dio flux modulation that it exhibits does not show typical 
nulling properties as displayed by some old isolated pul- 
sars; specifically, we have found strong evidence that its 
variation is highly frequency dependent. Although our 
orbital modulation analysis does not show a significant 
correlation between orbital phase and flux, the modula- 
tion could still be caused by a process related to a wind 
from its companion, which results in short time scale 
variations grouped in extreme modulation "events" . Ad- 
ditional monitoring of both the pulsar and of its com- 
panion may prove useful in this regard. 

We pointed out that this pulsar has an unusual combi- 
nation of characteristics: long spin period, very low-mass 
companion, high magnetic field and short orbital period, 
that are unexplained by standard binary pulsar evolu- 
tion scenarios. We propose that PSR J1744— 3922, along 
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with a several other binary pulsar systems, are part of a 
new class of low-mass binary pulsars which failed to be 
fully recyled. Specifically, we suggest three alternative 
scenarios for this class of binary pulsars. Distinguishing 
among them may be possible by improving our knowl- 
edge of the nature of their companions. We also reported 
the detection of a possible near-infrared counterpart to 
PSR J1744— 3922's companion, however, determining its 
nature will require detailed near- infrared/optical follow- 
up. 
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TABLE 1 

Orbital and timing parameters for PSR J1744— 3922 a 



Parameter Value 



Orbital Period, P orb (days) 0.19140635(1) 

Projected Semi-Major Axis, x (lt-s) 0.21228(5) 

Orbital Period Derivative, |P or6 | (ss" 1 ) <2x 10~ 10 

Projected Semi-Major Axis Derivative, |x| (lt-ss — ) < 7 X 10 — 12 



Derived Parameters 

Eccentricity, e <0.001 

Mass Function, /i (M Q ) 0.0002804(2) 

Minimum Companion Mass, M c (M ) > 0.085 

Surface Dipole Magnetic Field, B (G) 1.7 x 10 10 

Spin Down Energy Loss Rate (ergs -1 ) 1.2 X 10 31 

Characteristic Age, (yr) 17 X 10 9 



Note. — Numbers in parentheses represent twice the formal errors 
in the least significant digits as determined by TEMPO after scaling the 
TOA errors such that the reduced-x 2 of the fit was unity. The pulsar 

is assumed to have mas s 1.4 Mq. 

a Values reported from Ransom ct al] din prep.1) . 



TABLE 2 
Details of the radio observations 



^center 


BW b 


Num. chan. c Sampling T sys 


Gain 


rpe 


Num. obs. 


MHz 


MHz 


fis K 


K/Jy 


min 




Parkcs Telescope 


680/2900 g 


56/576 


256/192 500/500 68/31 


0.625/0.59 


20/20 


5 


1375 


288 


96 500 28 


0.71 


17 


69 


1400 


576 


192 250 32 


0.71 


15 


13 


Green Bank Telescope 


820 


48 


96 72 37 


2.0 


18 


4 


1400 


96 


96 72 24 


2.0 


25 


19 


1850 


96 


96 72 22 


1.9 


60 


3 


1950 


600 


768 81.92 22 


1.9 


210 


4 


4600 


800 


1024 81.92 20 


1.85 


257 


1 



Note. — Receiver temperatures and gains are estimated operating values. Parkes values are 
provided by J. Reynolds (private comm.). The system temperature corresponds to the sum of the 
receiver temperature and the sky temperature, which is determined from the 408 MHz all-sky survey 
and converted to other fre quencies by assuming a power-law spectrum having a spectral index —2.6 
l IHaslam et al.|[l982l. [1995T) . 

a Central frequency of the receiver. 13 Observing bandwidth. c Number of frequency channels. d System 
temperature. e Average total integration time per observation.' Number of observations^ Observa- 
tions were made simultaneously at these two frequencies. 



TABLE 3 

Percent of observations with no detection of 
PSR J1744-3922 



Frequency 


Parkcs 


GBT 


Both 




MHz 


% 




% 




% 




680 a 


100 


(5) 




(-) 


100 


(5) 


820 




(") 


75 


(4) 


75 


(4) 


1400 


33 


(72) 


32 


(19) 


33 


(91) 


1850 




(-) 





(3) 





(3) 


1950 




(") 





("1) 





(1) 


2900 a 


50 


('1) 




i 


50 


(4) 


4600 




(-) 





(i) 





(1) 


Total 


38 


(81) 


29 


(31) 


35 


(112) 



Note. — Numbers in parentheses represent the 
total number of observations for each band. 
a Observations at 680 MHz and 2900 MHz were si- 
multaneous. Excessive RFI contamination prevents 
us from using one of the 2900 MHz observations. 
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TABLE 4 

Estimated Pulsed Flux Density of PSR J1744— 3922 



Frequency 


Num. data a 


Average 


Standard Deviation 


Maximum 


MHz 




mjy 


mjy 


mjy 


680 


40 


<0.07 b 






820 


20 


0.10 


0.06 


0.22 


1400 


2244 


0.11 


0.16 


0.45 


1850 


142 


0.11 


0.04 


0.20 


1950 


852 


0.08 


0.04 


0.19 


2900 


40 


0.09 


0.08 


0.16 


4600 


258 


0.006 


0.012 


0.03 



Note. — Pulsed flux density values were derived using the radiometer 
equation, implicitly assuming that the off-pulse RMS is a good estimate of 
the system te mperature and that t he sk y emits according to the 408 MHz 
all-sky survey I IHaslam et al. 1982, 199|). Relat ive errors are estimated to 
be ~ 30 %. 

a Number of data points used at each frequency. The time resolution is 1 
minute per data point. b 3cr upper limit. 



TABLE 5 

Characteristics of partly recycled binary pulsars (P s > 8 ms) in 
the Galactic field in low-eccentricity orbits (e < 0.01) and 

HAVING LOW-MASS COMPANIONS (M Cimin < 0.2 M©) 



Name 




P 


logB 


Porb 


M • a 


Type 








ms 


G 


days 


M 






Case A 




<8 


8-9 




< 0.45 


He 


WD 


Case B 






9-10 




> 0.45 


CO 


WD 


PSR J1744- 


-3922 


172.44 


10.22 


0.19 


0.08 




? 


PSR B1718- 


-19 b 


1004.03 


12.11 


0.25 


0.12 


Bloated 


MS C 


PSR B1831- 


-00 


520.95 


10.87 


1.81 


0.06 




? 


PSR J1232- 


-6501 


88.28 


9.93 


1.86 


0.14 




? 


PSR J1614- 


-2318 


33.50 


9.14 


3.15 


0.08 




? 


PSR J1745- 


-0952 


19.37 


9.51 


4.94 


0.11 




? 


PSR B1800- 


-27 


334.41 


10.88 


406.78 


0.14 




? 


PSR J0407H 


-1607 


25.70 


9.16 


669.07 


0.19 




? 



a M c m i n refers to the minimum mass of the companion corresponding 
to an orbital inclination angle of 90° and assuming a mass for the pulsar 
of 1.3 5 M w . b In globular cluster NGC 6342. c [janssen fc van Kerkwiikl 
(2005) 
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Fig. 1. — Sample folded intensity profiles of PSR J1744— 3922 as a function of time for two 1400 MHz observations at Parkcs with 576 MHz 
bandwidth (panel A & B) and a 1950 MHz observation at GBT with 600 MHz bandwidth (panel C). The grayscale represents the intensity 
of the signal, with darker regions being brighter. The center of the gaussian-likc pulse profile should appear at the pulse phase = 0. 
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Fig. 2. — Sample light curves of the pulsed radio flux density as a function of orbital phase. Each data point represents 60 sec of data 
and orbital phases are defined so that 0.25 is when the companion is in front of the pulsar. Panel A, B and C are the two 1400 MHz 
observations and the 1950 MHz observation, respectively, displayed in Figure [T] We note that the flux drops below the detection limit at 
1400 MHz several times whereas it appears to be always above this threshold at 1950 MHz. Also, as illustrated in Figure [3] the variability 
is much stronger at lower frenquency. 
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Fig. 3. — Distribution of the measured pulsed flux density values at 1400 MHz (solid line) and at 1950 MHz (dashed line). The 3<r 
sensitivity threshold is ~ 0.02 and 0.01 mjy at 1400 and 1950 MHz, respectively. We observe that PSR J1744— 3922 rarely disappears at 
1950 MHz whereas there is a significant number of non-detections, centered about zero, at 1400 MHz (negative values are reported when 
the pulsed flux density is below the telescope sensitivity, meaning the flux determination algorithm has fit noise). Also, we note that the 
pulsed flux density is more variable and spans over higher values at 1400 MHz than at 1950 MHz. 
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Fig. 4. — Ratio of pulses detected with respect to the total number observed at different orbital phases for the 1400 MHz Parkes and GBT 
combined (panel A) and a similar plot showing the average pulsed flux intensity (panel B) for the Parkes 1400 MHz data (solid line) and 
the GBT 1950 MHz data (red dashed line). Orbital phases are defined so that 0.25 is the pulsar's superior conjunction (i.e. the companion 
passes in front of the pulsar). The scale of the error bars represents the estimated errors without scaling from the results of the Monte 
Carlo simulations, hence they are likely underestimated. 
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Fig. 5.— Near-infrared image of the field of PSR J1744-3922 in the K' band, obtained with AOB KIR at CFHT. The 0'.'34 positional 
error circle (3<x confidence) is shown, with the proposed counterpart at the centre. 
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Fig. 6. — Cooling tracks for He WDs made using the evolution sequences of [Dricbc ct aQ 1)19991 ) and our atmosphere models (see 
IBergeron eTaniT99al2l)0ll ; lBerg"eronll2001l . for more details). Lines show the cooling for constant masses of 0.179, 0.195, 0.234, 0.259, 0.300 
and 0.331 Mq, from bottom to top, respectively. The shaded region is the restricted range of absolute K' band magnitude inferred from 
the CFHT data, and the dotted vertical line indicates the characteristic age of PSR J1744— 3922. 
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Fig. 7. — Orbital period versus companion mass for binary pulsars in the Galactic field in circular orbits (e < 0.01). Symbols indicate the 
companion median mass (corresponding to i = 60°) and are coded according to the minimum mass: M c < 0.2 Mq, 0.2 Mq < M c < 0.35 Mq 
and M c > 0.35 Mq are black triangles, blue stars and red circles, respectively. Error bars cover the 90 %-probability mass range for randomly 
oriented orbits having i = 90° to i = 26°. The curves are the predicted P or (, — M c relationships for different metallicity progenitors (from 
Tauris fc Savoniielll999l) . PSR J1744— 3922 is identified by a square outline and pulsars listed in Table[5]with triangle outlines. The plot 
also includes the globular cluster pulsar PSR B1718— 19 (marked by a circle outline) because it resembles PSR J1744— 3922 (see Section 
12.41 1. Data from the ATNF Pulsar Catalogue (Manchester ct al. 2005) (http://www.atnf.csiro.au/research/pulsar/psrcat/). 
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Fig. 8. — Inferred surface dipolar magnetic field strength versus spin period (panel A) and spin period versus orbital period (panel 
B) for binary pulsars in the Galactic field in circular orbits (e < 0.01). Symbols are coded according to the companion minimum mass: 
AI C < 0.2 Mq, 0.2 Mq < M c < 0.35 M Q and M c > 0.35 M are black triangles, blue stars and red circles, respectively. PSR J1744-3922 is 
identified by a square outline and pulsars listed in Table [5] with triangle outlines. The shaded area is the approximate region where is the 
proposed class of binary pulsars similar to PSR J1744— 3922. The plots also include the globular cluster pulsar PSR B1718— 19 (marked 
by a circle outline) because it resembles PSR J1744— 3922 (see Section 12.41 1. The dashed line is the best-fit for a power-l aw, B oc P a , with 
a = 1 .13. We excluded PSR B1718— 19 from the fit as it is in a globular cluster. Data from the ATNF Pulsar Catalogue (Mancheste r et aTl 
2005) (http:/ /www. atnf.csiro.au/research/pulsar/psrcat/ ). 



